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Abstract 

In this paper we present a study of an exciton system where electrons and holes are confined in 
double quantum well structures. The dominating interaction between excitons in such systems is 
a dipole - dipole repulsion. We show that the tail of this interaction leads to a strong correlation 
between excitons and substantially affects the behavior of the system. Making use of qualitative 
arguments and estimates we develop a picture of the exciton - exciton correlations in the whole 
region of temperature and concentration where excitons exist. It appears that at low concentration 
degeneracy of the excitons is accompanied with strong multi-particle correlation so that the system 
cannot be considered as a gas. At high concentration the repulsion suppresses the quantum degen- 
eracy down to temperatures that could be much lower than in a Bose gas with contact interaction. 
We calculate the blue shift of the exciton luminescence line which is a sensitive tool to observe the 
exciton - exciton correlations. 
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FIG. 1: Excitons in coupled quantum wells, d is the separation between the centers of the wells 



(given d the dependence of the exciton - exciton interaction on the well widths is very weak 16l|). 
d and the average separation between excitons are assumed to be larger than the exciton radius 
ax- 



INTRODUCTION 



A very active investigation of excitons in coupled quantum wells for more than two decades 
was first motivated by the possibility to reach Bose condensation and superfluidity in this 
system. Further experiments discovered a very large number of related phenomena and quite 
rich physics of the system (see Refs. and references therein)- Theory predicts the 

existence of many phases with different and unusual properties. 

mm 

The most expected 

and searched for is the settling in of coherence in such 2D exciton systems. A coherence of 
the exciton Bose condensate has to reveal itself in some coherent properties of the exciton 
uminescence. Investigation of the luminescence led to discovery of not only its coherence 
8, 9] but also a number of patterns not completely understood so far 
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,ll2,ll3|. 



A substantial role in these phenomena is played by the interaction between excitons. 
Typically, interaction between bosons (^He atoms and alkali-atoms) is of a short range and 
the theory of non-ideal Bose gas has been developed for contact interaction, jl^ . Iisl In 
coupled quantum wells where the electrons and holes are separated in the two adjacent 
layers, all the indirect, dipolar excitons that are formed by the coulomb binding of pairs 
of these spatially separated electrons and holes are polarized in the same way and their 
interaction is mainly dipole - dipole repulsion, FiglH Contrary to the contact interaction 
the dipole - dipole interaction has a signiflcant tail and due to this tail the exciton gas in 
some respects is dramatically different from Bose gas with contact-like interactions. 

The dipole - dipole repulsion increases the exciton energy and leads to a blue shift of 
its luminescence line. The blue shift is usually evaluated as an average value of the exciton 
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interaction energy with other excitons and can be easily obtained with the help of the mean 



field approximation that produces "plate capacitor formula" 



17| 



„ Anne'^d 

Eint = , 1.1 

K 

where n is the exciton 2D concentration, d is the separation between the centers of the 
wells, and k is the dielectric constant. This formula can be understood in the following 
way. Equal concentration n of electrons and holes in two wells creates potential difference 
A(f) = ATTned/n between them. Creation of one more indirect exciton in these wells requires 
transfer of an electron or a hole from one well to the other. In the presence of other excitons 
this increases the necessary energy by eA0 that gives Eq. fll.ll) . This expression is typically 
used in experiments for an estimate of the exciton concentration from the measured blue 
shift of the luminescence. 

Recently Zimmermann and Schindler ISj noticed, however, that dipole - dipole repulsion 
leads to a significant exciton pair correlation. The repulsion prevents excitons to come very 
close to each other and creates a depleted region around each exciton. The pair correlation 
leads to a reduction of the coefficient in Eq. fll.ll) by about 10 times (depending on the 
excitons temperature). [iS^ This means that previous experimental estimates of the exciton 
concentration based on the capacitor formula underestimated the concentration by up to 10 
times! 

It makes sense to note that the reduction of the pair correlation function to zero at small 
distance in a Bose gas with repulsion has been noticed long ago and used in the construction 



of a many particle variational wave function 
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20| It is well known in exact solutions for 



one dimensional Bose gas. |2lL 1221 . l23l . l24l | In the exciton gas with dipole - dipole repulsion 



this behavior of the pair correlation function was noticed by Astrakharchik et al. 71]. This 
behavior was also used by Kash et al to explain a narrowing of the exciton luminescence 
line at low temperature, jisl 

The significant pair correlation between excitons substantially affects our understanding 
of the behavior of this system. This motivated us to develop a general picture of the interac- 
tion and correlation between excitons in the whole temperature - concentration plane. This 
is the subject of this paper. We show that Eq.f ll.ll) has a very limited region of applicability 
and limited accuracy. With a decrease of the temperature and an increase of the concentra- 
tion the role of the pair correlation becomes crucial for interpretation of such phenomena 
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as the blue shift. It appears also that the degeneracy of the exciton system at any small 
concentration is accompanied by setting in of not pair but multi-particle correlation charac- 
teristic for liquids. In other words, a consistent theory of a quantum coherent state has to 
include multi-exciton correlation at any dilute gas. More than this, at higher concentrations 
the dipole - dipole repulsion reduces the overlap of the wave functions of different excitons. 
As a result the exciton system can be considered as classical liquid down to temperature well 
below than the temperature of quantum degeneracy of a system with contact interaction. 

Our main purpose is to develop a qualitative understanding of the structure of the system 
of interacting bosons. Therefore we don't pursue a high accuracy of the results but instead 
use simplified models and qualitative arguments. Although the results of such an approach 
are really accurate only in extreme cases with respect to some large or small parameters, 
they allow us to make analytical calculations and produce a very clear physical picture of 
relevant phenomena. Availability of such a picture facilitates precise calculations when they 
are necessary. 

In the next section we consider in detail the applicability of the mean field approximation 
that neglects any correlation between excitons. An exciton gas at lower temperatures where 
strong pair correlation is important but quantum correlations are still negligible is considered 
in Sec JIIII Further reduction of temperature at low exciton concentrations, when the exciton 
wavelength becomes larger than the characteristic scale of exciton - exciton interaction, leads 
to an important role of quantization in the exciton - exciton scattering while the exciton gas 
itself is still statistically non-degenerate. This situation is considered in SeclIVI In SecJVl 
we consider the situation where multi-exciton correlation is important. This happens at 
low enough temperature in a dilute system and in a more wide range of temperatures in a 
dense system. It appears that with a temperature decrease in a dilute system degeneracy 
is accompanied by a set in of a multi-exciton correlation. Contrary, in a dense system the 
classical multi-exciton correlations appear to be more important for the blue shift than the 
quantization of exciton dynamics. In the last section we discuss the obtained results. 

The discussion of exciton correlation in different parts of the concentration - temperature 
plane involves quite a large number of physical parameters. To facilitate reading of the 
paper we present here the list of main parameters with their physical definitions. Exact 
mathematical definitions, if necessary, will be given as soon as the corresponding parameter 
comes into the discussion. 
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TABLE I: The main parameters characterizing an exciton gas in coupled quantum wells 



d 


Separation between electron and hole wells 


CLB 


Bohr radius with a reduced electron - hole mass 


nip, mh 


Electron and hole masses 


b 


Bohr radius with an exciton mass 


m 


Reduced electron - hole mass 


ax 


Exciton radius 


M 


Exciton mass 


ro 


Average minimal distance between excitons 




Exciton binding energy 


kT 


Thermal wave vector of an exciton 


n 


Exciton concentration 


U{r) 


Interaction energy between excitons 


Eint 


Average interaction energy of a single exciton 


9{r) 


Exciton pair correlation function 



II. MEAN FIELD APPROXIMATION 



The interaction energy between excitons in coupled quantum wells is 



f/(r) 



2e^ 



1 



Vr2 + 6/2 



(2.1) 



According to Refs.j4, 18| the attractive Van der Waals interaction and the exchange in- 
teraction are small in practically important values of d. The simplest way to obtain the 
average interaction energy is to assume that all excitons are distributed randomly and in- 
dependently of each other with an average concentration n, see Figj2l Then the average 



number of excitons in an area element d'^r is nd'^r and the average interaction energy is 



E.. 



int 



U {r)nd'^r 



Aime^d 



(2.2) 



This result means that the main contribution to Eint comes from the interaction between 
excitons at distance of the order of d: {e^ / K.d){mtd^) = ime^d/K,. 

Eq. (]2.2p is valid only under a few assumptions. The first is that the average distance 
between excitons has to be larger than d. In other words, the concentration cannot be too 
large, 

nrf^ < 1 . (2.3) 

In the opposite case the electron - electron and hole - hole repulsion is stronger than the 
electron - hole attraction and it is hardly possible to expect a stable exciton phase of the 



system. [26 
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28|, 



29| In all following considerations we assume that inequality Eq. (12.31) is 



met. Practically, this limitation in coupled well structures with (i ^ 10 nm means than the 
concentration has to be much smaller than 10^^ cm"^. 
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FIG. 2: Uncorrelated motion of excitons. The radius of each circle depicts the characteristic scale 
of the repulsive potential d. A random motion takes place only at low concentrations, Eq. (|2.3p . 
and rather high temperatures, Eq. (|2.5|) . 

The other main assumption used in the derivation of Eq. fl2.2l) is the absence of any 
correlation between excitons. In reality, if the kinetic energy of the relative motion of any 
two excitons is E they can approach each other only to a distance larger than ro{E) where 
Tq is the root of the equation 

f/(ro) = E . (2.4) 

The correlation between excitons can be neglected only if the region where the correlation is 
important is very small, i.e., tq <C d. If the temperature of the exciton gas is T then E T 
and this condition can be written as 

T > e^/nd . (2.5) 

Condition (12.51) is equivalent to kx ^ \/\/bd where = y/2MT/fi is the exciton thermal 
wave vector and b = h^n/Me^. The expression for h differs from the Bohr radius as = 
h^K/me^ only by replacement of the electron - hole reduced mass m = niemh/ {nie + mh) 
with the exciton mass M = rrie + rrih. In GaAs/AlGaAs structures the electron effective 
mass rrie = 0.067 and the hole effective mass at the bottom of hhl subband in a quantum 
well rrih = 0.14. This gives m ^ 0.045, M ^ 0.21, and therefore ~ 14 nm, and 6 ~ 3 
nm which allows us to assume in further calculations that h <^ d. Therefore the condition 
(12. 5p means also that kx ^ l/(i, i.e., the exciton wavelength is much smaller than the 
characteristic length scale of the potential. This justifies a classical consideration of the 
interaction between excitons. Also, this inequality in combination with Eq. (l2.3p leads to the 
inequality k^ ^ n which means that the exciton thermal wavelength is much smaller than 
the average inter-particle distance and therefore the exciton gas is non-degenerate. 
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t makes sense to note that a quantum mean field approximation also leads to Eq. (12. 2 p . ly, 
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32j However, if the exciton radius ax is of the order of or larger than d the exchange 



interaction also appears to be important. 



161 ] Variational calculations for GaAs/AlGaAs 



structures with an exciton wave function ip{r) = Ae"^"'^+'"^/^"^ yields, for layers separation 
of d = 10 nm, 12 nm, and 14 nm the corresponding exciton radii of ax = 8.7 nm, 9.3 nm, 
9.7 nm, and exciton binding energies of 4.7 meV, 4.2 meV, and 3.8 meV respectively. That 
is although d > ax this inequality is not very strong and there are quantum corrections to 



the interaction in fl2.1D. 
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The exciton binding energy also puts an upper limit to the temperature where the above 
mean field description is viable because at temperatures of the order of or larger than the 
binding energy the majority of the excitons dissociate. According to Eq. (12.51) the low limit 
for the temperature for d = 10 nm, 12 nm and 14 nm is 54 K (4.7 meV), 48 K (3.8 meV) 
and 44 K (8.6 meV) respectively. 

The bottom line of these estimates is that Eq. ( 12. 2p is a rather poor estimate: it is really 
valid only in a the temperature range where a significant part of excitons is dissociated. 

To conclude this section we comment on a usage of the mean field approximation in 
explaining phenomena other than the exciton luminescence blue shift. The interaction of 
one of the particles with all others can be described with the field created at the particle by 
the environment. In general, this field fluctuates in time and from particle to particle due to 
different dynamics of particles creating it. If the particle interacts simultaneously with many 
others and they are not correlated then these fluctuations are cancelled and their resulting 
amplitude is much smaller than the average value of the field. This is the foundation that 
makes the mean field approximation valid. In the exciton gas with dipole - dipole interaction 
the situation is quite different. The average interaction energy Eq. (l2.2p is much larger than 
the interaction between excitons at average distance between them: f/(n'^/^) ~ e^ip'rc'l'^ j k 
and U{n^^^'^) / Eint ~ {nd'^Y^'^ <^ 1. This means that the main contribution to Eint comes 
from rear pairs of excitons with the distance much smaller than the average one. The large 
amplitude of the field fluctuations is also confirmed by the calculation of the average of the 
interaction energy squared: 



f/2 = y U\r)nd^r = 2tx In y • (2.6) 

Discarding the logarithm that comes from the cutoff of the minimal distance between exci- 



tons we see that El^^/V^ ~ nd^ <^ 1. (Under this condition characterizes the lumines- 
cence hne width, but a discussion of this point is beyond the scope of the paper.) 

The message is that while the mean field approximation for the average interaction gives 
a correct result in the mentioned range of parameters, the calculation of other quantities in 
this approximation can lead to large errors. It is also important to keep in mind that in 
spite of the purely classical arguments this statement is true also in the quantum limit. 

III. CLASSICAL EXCITON GAS 

When the temperature goes down at some point it becomes smaller than the Coulomb in- 
teraction at distance d that violates the condition of the mean field approximation, Eq. (]2.5p . 
and the correlation of excitons cannot be neglected anymore, Figj3l At T < / nd the rele- 
vant parameter characterizing the scale of the interaction potential is not d but the minimal 
distance between excitons tq because tq > d. At the temperature where tq crosses the value 
of (i, condition (12.31) leads to 

nrl < 1. (3.1) 

This means that right below the temperature T ~ e^/ztd there exists a region where the 
average distance between excitons is larger than ro and it is possible to take into account 
only pair correlations because the probability to find three excitons in mutual proximity is 
negligible. The condition tq ^ d (i.e., T ~ / nd) means also that /ct'^o ~ k^d ^ kxVdb ~ 1 
(the last two relations follow from d ^ b and T ~ e'^/nd respectively) and the interaction 
between excitons can be still considered classically. The same relations mean that ^ n, 
i.e., the exciton gas is non-degenerate. A growth of tq compared to the exciton radius also 
makes the exchange corrections to the interaction (12.11) less important. As a result at the 
temperature region right below e^/nd it is possible to consider excitons as classical particles. 
The blue shift in this temperature and concentration region can be evaluated as 

Eint = nj U{r)g{r)d^r , (3.2) 

where g{r) is the pair correlation function. For any given exciton, ng{r)d'^r is the average 
number of excitons within area d'^r at distance r from it. According to this definition 
g{r)\r^oo = 1 because at large distance any correlation between excitons disappears. In the 
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FIG. 3: At low temperatures, T < / Kd, the exciton kinetic energy cannot overcome the repulsion 
between the excitons. The minimal distance between excitons in average is larger than ro defined 
in Eq. (j2.4p where the energy E ^ T. The correlations between excitons cannot be neglected. 



leading order in htq 



9[r) 



-Uir)/T 



see, e.g., Ref.|33| Sec. 32]. Substitution of Eg. (13. 3 p into Eg. (13.21) gives 



int 



(3.3) 



(3.4) 



With the interaction Eg. (l2.ip this expression is reduced to a function of only one parameter. 



E., 



int 



fE 



Function fsix) is plotted in FigJH The asymptotical behavior of this function is 



(3.5) 



1 + - In - + C - 1 

X \ X 



F(4/3) , 



X 



/3 



X > 1 , 

x -C 1 , 



(3.6) 



where C ~ 0.577 is the Euler constant. At high temperature, Eg. (l2.5p . the expression for 
Eint (13. 5p is reduced to the mean field expression of Eg. (12.21) . In the opposite limit, where 
ro ^ d, the interaction potential is simplified: 

e^d^ 



U{r) 



r:^d 



(3.7) 



and Eg. (l3.4p gives 



'e2^2\ 2/3 

Eint = 27rF(4/3)n ( T^'^ = 2TiV{A/?>)nrlT 



(3. 




FIG. 4: A plot of dimensionless function fsix) of Eq. (j3.5p that defines the blue shift at low 
exciton concentrations, Eq. ()3.ip . and moderate temperatures, Eq. ()3.14p . 

where tq = (e^cP /kTY^^. Qualitatively this result can be understood in the following way. 
Around each exciton there is a depletion region with a radius ~ tq. Without repulsion this 
region would contain rarg excitons with average energy T. The energy necessary to force all 
of them out of the region is of the order of nr^T. 

One can expect that is of the order of the first virial correction to the chemical 
potential of the exciton gas. This can be easily checked. The first two terms of the virial 
expansion give (Ref. |33l] Sec. 23) 



n 



2mT 



/2mT 



r . 



Solution of this equation with respect to the chemical potential C leads to 



C = Co + AC 



Co = Tin 



and 



AC = -nT J [e-^(-^/^ - l] d'r 



2MT 
2Tine^d 



/c 



exp 



2 /I 1 



KdT 
tdt 



In extreme cases 



AC 



^r(l/3)(^l nT^'^ 



T > 



nd 

T < — . 

Kd 



(3.9) 

(3.10) 

(3.11a) 
(3.11b) 



(3.12) 
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That is Eint and have the same dependence on parameters and differ only by a constant 
factor. This difference results from the difference of the definitions: Eint is the correction to 
the average energy released in an exciton recombination while A( is the correction to the 
energy released when an exciton is removed without a violation of the equilibrium in the 
exciton gas. 

Finally, we note that in the case of T <^ e'^/nd when Eint is given by Eq. fl3.8p . the 
interaction at the average distance U{n~^^'^) is still small compared to Eint- U{n~^^'^) / Eint ~ 
n^/^ro ^ 1. The mean field approximation then gives correct expression only for the average 
energy because 



i.e., the fluctuations of the interaction energy are larger than its average value. 

The results of this section, Eqs. fl3.4p and (13.81) . are valid under two conditions: that of 
small concentration, Eq. (l3.1l) . and 



that validates classical description of the interaction. 

The parameter utq is the gas parameter which is the ratio of the exciton interaction 
energy to its kinetic energy. The same parameter indicates the strength of the exciton - 
exciton scattering. If the impact parameter in a scattering event of two excitons is ^ tq then 
the scattering angle is large. Therefore the scattering cross-section (in 2D case it has units 
of length) is ~ tq. Respectively, the mean free path is / ~ l/nrg. The scattering is weak, 
i.e., three and more particle scattering can be neglected if / is larger than the interparticle 
distance which means that the gas parameter is small, Eq. (l3.1l) . If this condition is violated 
the exciton system cannot be considered as a gas, it is a liquid. On the other hand, parameter 
kxTQ characterizes quantum corrections to scattering of excitons. Eq. (l3.14l) is stronger than 
the non-degeneracy condition because fcr'^o = {kx / n^^'^){n^^'^rQ) <^ k^lvt^^'^ because the 
characteristic scale of the potential tq is smaller than the average distance between excitons. 



When temperature goes down further both condition of small concentration, Eqs. (l3.1l) . 
and classical description of exciton - exciton scattering, Eq. (l3.14p . are at some point violated. 
Which one is violated first depends on the concentration. If n < {h/2(PY then quantum 




(3.13) 




(3.14) 



Eq.iO). 
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effects in the exciton scattering become important when the exciton system can still be 
considered as a non-degenerate gas. In the opposite case with reduction of the temperature 
the exciton system becomes a liquid before any quantum corrections to the scattering process 
become pronounced. 

IV. EXCITON GAS WITH QUANTUM SCATTERING 
If 

n < {h/2d^f (4.1) 

and kTTQ ^ 1 which is equivalent to 

kT < h/2d\ (4.2) 

then quantum corrections to the exciton - exciton scattering are important but the exciton 
gas is yet non-degenerate until kx ^ n^/^. For the calculation of the interaction energy it is 
then still possible to use Eq. fl3.2p but g{r) has to be modified to include quantum corrections. 
This can be done in the following way. 

A wave function describing a state of two excitons can be factorized into a wave function 
of the center of mass and a wave function ip{r) describing their relative motion (see Appendix 
'il){r) is characterized by a few quantum numbers but in equilibrium the occupation of a 
state depends only on its energy. This means that, given the energy of the relative motion 
E of any two excitons, the probability to find one exciton at the distance r from the other 
is {\i'{'r)\^) E where (. . . )_e is the average over all quantum numbers (e.g., the direction of 
the wave vector) except E. For a non-degenerate exciton gas the probability density for an 
exciton to have energy E is {1/T)e~^/'^ . That is 

1 

g{r) = - {mr)\')Ee-^/^dE (4.3) 

where ipi''^) has to be normalized in such a way that g{r)\r^oo = 1 which corresponds to 
usual normalization for scattering problem. 

In a semiclassical approximation, when the exciton wavelength is smaller than the length 
scale of the interaction potential, krro ^ 1, 

^<-' = W^WW ' 
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FIG. 5: A sketch of an exciton gas with quantum scattering. The radius of the white circle is of 
the order of the exciton wavelength, ~ l/kx, the radius of the gray circle is the classical minimal 
distance between excitons, ~ tq, and the radius of the black circle is the quantum minimal distance 
between excitons, ~ d? /b. Here the exciton wavelength is larger than ro and the scattering is 
essentially quantum mechanical. On the other hand the average distance between excitons is 
larger than their wavelength so that the gas is non-degenerate. 

where S satisfies the equation 

^ = E- U{r) . (4.5) 
Substitution of Eq. fl4.4p in Eq.f l4.3p results in 



which is identical to Eq. (13.31) for = ^JtJtx. Fig. 11 of Ref. 18] for g{r) is related to an 
intermediate case where ^t^^o ~ 1-7 and it differs from our classical expression by 20% in 
the scale of r which comes from not very large value of fcr'^o and our simplification of the 
interaction between excitons. 

If /ct^o ^ 1 the result strongly differs from the classical case. The wave function V'(^) 
penetrates under the repulsion barrier and the minimal distance between excitons is char- 
acterized not by ro anymore but rather by the distance at which '?/'(r) falls off. The dipole 
- dipole repulsion e^c?^/ kt^ is trying to push the wave function to larger r while the kinetic 
energy ~ {h^ /Mr'^) is trying to spread it to all available space and in particular to as small 
values of r as possible. The distance at which ip{r) falls off is characterized by the same 
order of magnitude of these two tendencies. This gives the distance of ~ d? /h, FigJSl 

When the interaction potential is approximated by interaction between point like dipoles, 



Eq. (l3.7l) . this distance is smaller than tq: ro/{d'^/b) ~ {b/krd'^y^^ > 1 because of Eq. (14.21) . 
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But due to 6 ^ (i it is much larger than d which justifies the approximation of point-hke 
dipoles, Eq. fl3Tl) . 

The interaction energy can be estimated as U{d'^/b) (under the condition kxTQ <^ 1 it 
is larger than T) multiplied by the number of excitons in the important region, n{d'^/b)'^, 
that gives Eint ~ h'^n/m. The exact calculation (Appendix [XI) shows that ipi^) off 
exponentially when r < d'^/b and is a relatively slow function of r when r > d'^/b. On 
the other hand, the potential U{r) falls off as when r > d'^/b. As a result, the main 
contribution to the interaction energy Eg. (13.21) comes from the region r ~ d^/b. In this 
region the wave function is given by Eqs.( lA7l) and ( ]A12I) : 



where k = yME/h. The logarithmic dependence of this wave function on the energy is 
very weak and with a logarithmic accuracy the substitution of Eq. (l4.7p into Eq. (l4.3p leads 
to the following expression for the correlation function: 



As a result, 



87rn e^d^ . f 2d \ dr 27rh^n , , 



hi^{kTdM>) K Jq " \ Vb^J Mln^krdyb) 



where Eq. (6.576.4) of Ref. |35l| has been used. Note that in this case again Emt is 



larger than the interaction at the average distance between excitons: U{n ^/'^)/Eint ~ 
{d'^n^/yb)\n^{kTdyb) < 1 which is due to Eq.dH]). 

The gas parameter in the quantum case is different from classical one, nrg. Exciton - 
exciton scattering is strong if two excitons approach each other to a distance equal to their 
wavelength. In other words, the scattering cross-section is of the order of the wavelength with 
accuracy of a logarithmic correction (this is the well known difference between the 2D and 



the 3D case where at small wave vectors the cross-section goes to a constant, see Ref. |34| . 
Problem 7 to Sec. 132, and Appendix in spite of the fact that the length scale of the 
potential d'^/b is smaller than the wavelength. Respectively the mean free path of excitons 
is / ~ krp/n. The gas condition in quantum case is the absence of correlations between 
different scattering events which means that the wavelength has to be much smaller than 
the mean free path, i.e., kTl ^ 1 or 

n <^ k^ln\kdyb) . (4.10) 
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FIG. 6: A rough sketch of quantum exciton hquid. The radius of the white circle is of the 
order of the exciton wavelength, ~ I/^t, the radius of the gray circle is the classical minimal 
distance between excitons, ~ ro, the radius of the black circle is the quantum minimal distance 
between excitons, ~ The average distance between excitons is of the order or smaller than 

their wavelength and the system is degenerate. But the same condition means that multi-particle 
correlation is strong and this is a liquid. The picture cannot demonstrate that different excitons 
are not distinguishable due to an overlap of their wave functions. 

This inequality is identical with Eint/T <^ 1 and also with the condition of non-degeneracy 
with an accuracy of the logarithmic correction. Practically the logarithm is not very large 
and in the gas state the exciton system is non- degenerate while the degeneracy is accompa- 
nied with strong interactions and multi-particle correlations between the excitons which is 
characteristic for liquids. 

V. LIQUID STATE OF THE EXCITON SYSTEM 

In this section we consider the temperatures and/or concentrations beyond the limits 
specified in previous sections. In those cases multiexciton correlations are important and 
the problem is not reduced to a two-particle problem, Figj6l Actually all close neighbors are 
correlated although a long range correlation may not exist. This situation is characteristic 
for liquids and therefore we use the term "liquid" for such states of exciton system. 

An exact microscopical theory of liquids does not exist and to avoid poorly controllable 
and rather complicated approximations we make estimates based on reasonable physical 
arguments. These estimates give not only an approximate value of the blue shift but also 
its dependence on the concentration and temperature. Also, the arguments and estimates 
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allow us to develop a general understanding of the structure of the exciton system at low 
temperatures and high concentrations. 

The picture is different depending on the exciton concentration compared to (6/2(i^)^. 



A. Low concentration 

If the concentration is low according to quantum criterion Eq.f l4.ip . and the temperature 
is low enough so that multi-particle correlations take place, i.e., the inequality in Eq. fl4.10l) is 
violated, then Eq. (14.91) can be considered as a good estimate that can have only logarithmic 
corrections. Indeed, the temperature in Eq. (14.91) enters only in the argument of the logarithm 
and only this argument can change when the temperature goes down. We emphasize once 
again that violation of condition (14.101) leads not only to degeneracy but also to multi-particle 
correlation which makes the theory of dilute Bose gas unapplicable. 

We would like to attract an attention to a generally known fact that in a dilute 2D Bose 
gas the characteristic energy at low temperatures does not depend on the coupling constant, 
except logarithmic corrections. This comes from the virial theorem, i.e., from comparison of 
the interaction energy and kinetic energy and is a generalization of Eq. (l4.9p to any interaction 
between particles. In a dilute gas only two particles can be at the distance where their 
interaction is important (large quantum uncertainty of the distance ~ l/fcy means only that 
we cannot be sure that they are at this distance). If the interaction between particles is U{r) 
then from the virial theorem it follows that U{r) ~ h'^/Mr'^. The value of r obtained from 
this relation is the distance at which the interaction is important. The interaction energy 
of a particle is the interaction energy between two particles at the distance ^ r times the 
probability that two particles come to this distance, nr^. This gives nfi^/M. In a liquid 
the estimate can have a numerical factor characterizing the number of particles within the 
interaction radius. This energy gives a temperature scale for both Bose condensation 



and Kosterlitz - Thouless transition 



B. High concentration 



36 



38|. 



37| 



If the concentration is beyond its quantum limitation, i.e., Eq. (l4.ip is not satisfied then 
decrease of the temperature or increase of the concentration leads to violation of the classical 
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gas condition, Eg. (13. ID . while the system is still non-degenerate, i.e., Eq. (l3.14p holds. In 
this case the exciton system becomes a classical liquid. As long as Eq. (13.71) is valid the 
dimensional analysis allows us to express Emt a function of only one parameter: 

E^nt = ^ n^'f ( ^ n^/^ . (5.1) 



K \ kT 

According to Eq.([3l]) f{x) = 27ir{A/3)x-^/^ at x < 1. 

When [e^cp / kT)ii?^'^ = r^n^^"^ grows and becomes of the order of unity a free motion of 
excitons between collisions becomes impossible because each of them is confined in between 
its neighbors. In other words an exciton is in a highly excited state in a potential well 
formed by its neighbors. The size of the well is i? ~ n'^^"^ and this semi-classical picture is 
valid as long as the size of the confinement region is much larger than the exciton thermal 
wavelength, i.e. krpR ^ 1. The energy at the bottom of the potential well is ~ ze^d^ / {z is 
the number of nearest neighbors) and is of the same order as the depth of the well. Potential 
wells for different excitons are different, they are not static and sometimes some excitons 
overcome or tunnel across the surrounding barrier. But at utq > 1 these rear occasions do 
not affect the estimates. In general, this picture is similar to a simple classical liquid and 
the formation of the potential wells is the starting point of the formation of a short range 



order characteristic for liquids. 40|] We emphasize that we mean a formation of a short range 



order typical in liquids but not crystallization and a formation of a long range order [4l|], 
FiglZl Further reduction of the temperature brings particles to lower levels in the potential 
wells and makes the wells more stable. A stronger confinement of the wave functions of each 
exciton reduces their overlap. When the size of the exciton wave function becomes smaller 
than R, the potential for each exciton can be approximated as 

UH,ir)^'^(c^ + C,Qi, (5.2) 



where r is the distance from the minimum of the well. We estimate the constants Ci and 
C2 assuming a short range order, i.e., there is a crystal structure around an exciton within 
/ coordinate circles but beyond this region the exciton positions are not correlated. This 
assumption gives 
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FIG. 7: A rough sketch of a short range order in a classical liquid. In this example the order 
extends to two coordinate circles, and at larger distance excitons are not correlated. (Note that 
this picture is oversimplified: similar order exists around most of the excitons). 



TABLE II: Constants which characterize the effective potential well for each exciton in a liquid, 
Eq. (j5.2p . / = corresponds to very small correlations between the exciton positions, and / = oo 
corresponds to a crystal structure. 



Coordinate circle 


Square lattice, z = 4 


Hexagonal lattice, z = 6 


I 


Ci 




Ci 







6.28 


4.71 


7.25 


5.44 


1 


8.44 


10.67 


10.19 


14.55 


2 


8.55 


11.18 


10.78 


15.07 


3 


8.72 


11.29 


10.65 


15.08 


oo 


8.84 


11.37 


10.84 


15.13 



Here R = [nsin(27r/z)]~^/^ is the lattice constant which is the same as the radius of the 
first coordinate circle, Rj is the radius of the jth coordinate circle, and Zj is the number of 
particles at the jth coordinate circle, zi = z. The resulting values of the constants appear 
to be very weakly sensitive to the radius of the order, see Table [ITl (compare Ref. 42|). 

If T ~ [e^d? / k)!!?/"^ (i.e., nrg ~ 1) then r R and there is no short range order in the 
system. However, if T <^ [e^d? / K)rc'/'^ then the short range order does exist, and most of 
the excitons are at the ground state of their corresponding potential Eg. (15.21) . The energy 
of the ground state above the bottom of the potential well and its radius are 



2^2 



r e^d? 2C2 n'n d^n^l^ 



Tgs = R 



1 bR 



1/4 



(5.4) 



Mi?2 M \ b ' " \2C2 

The characteristic size of the exciton wave function r is controlled by the temperature or 
the energy of the ground state, whichever is larger, and in any case when T <^ [e^d^ / K)rr'/'^ 
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the size of a single exciton wave function is much smaller than the distance between the 
excitons. 

The inequality r <^ R allows us to make two conclusions. First, the bottom of potential 
Uiiq{r), Eq. fl5.2p . gives a good estimate for the interaction energy 

E^nt ~ . 5.5 

K 

Comparison of this expression with Eq.f l5.ll) shows that f[x-^ oo) ~ 10. This estimate 
does not include possible logarithmic corrections. 

Second, it is possible to estimate the overlap of the wave functions of adjacent excitons. If 
d = 12 nm and n = 2x 10^^ cm~^ the estimate according to the wave function in the harmonic 
potential of Eq. fl5.2p gives for the overlap a value of ~ 0.14. The actual value is even smaller 
because at r ~ i? the potential barrier is steeper than the harmonic one. Due to the small 
wave function overlap the temperature at which the phase or/and spin coherence ^ in the 
exciton system is set in is reduced compared to its the expected value ~ U^n/M. This points 
to a possible non-monotonic dependence of the quantum coherence onset temperature on 
the concentration, and it suggests that a lower density exciton system may become quantum 
coherent at higher temperatures than a higher density system, which is a-priori non intuitive. 
In other words, a long range interaction suppresses quantum degeneracy. 

At low temperature the physics of the transition between quantum liquid and classical 
liquid with growth of the exciton concentration is the following. At low concentration 
n ^ {b/d^y, Eq. fl4.ip . according to the exciton wave function Eq. fl4.7p around each exciton 
there is a circle with radius ~ c/^/fe inside which the wave function of any other exciton is 
exponentially small. However, the radius of this circle is much smaller than the average 
distance between excitons, n~^/^, and any exciton wave function can easily spread at the 
area that contains many other excitons avoiding their depleted circular regions. That is 
wave functions of different excitons overlap forming a quantum liquid (not a gas because 



of strong exciton - exciton scattering). 44| With growth of the concentration the average 
distance between excitons decreases that makes spreading of each exciton wave function to a 
wide area more difficult. This reduces the overlap of the wave functions of different excitons. 
Finally, when the average distance becomes smaller than the radius of the circular depleted 
region, -^/^ < b/d"^, the wave function of nearly each exciton appears to be confined in 
between its nearest neighbors. The overlap of the wave functions of adjacent excitons is very 
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weak and the system becomes a classical liquid. 



VI. DISCUSSION 

Estimates made in the previous sections open the possibility to develop a general picture 
that demonstrates a role of correlations in the exciton system at the whole n — T plane. 
This picture is presented in Fig. |H1 Correlations are not important and the mean field 
approximation is applicable only in region I. In region II the exciton system can be considered 
as a classical gas with strong pair correlations. In region III, contrary to region II, the exciton 
- exciton scattering is described by quantum mechanics. In other respects this region is 
similar to region II. Reduction of the temperature from region III to region IV leads to 
degeneracy of the exciton system. But simultaneously a strong multi - particle correlation 
is set up. The system cannot be considered as a dilute gas, and the mean free path does 
not exist. Rather surprising is the existence of region V where the system behaves as a 



classical liquid down to temperatures well below h^n/M (compare Ref. 43l]). The reason is 
that strong repulsion between excitons squeezes the wave function of each exciton to an area 
smaller than the average average area per one exciton. In this region a short range order 
appears and with further reduction of temperature its correlation radius grows. However, 
contrary to regular classical liquids, the attractive part of the exciton - exciton interaction 



is negligible 4, llSl] and it is likely that a long range order is settled not as a result of a phase 
transition but as gradual growth of the correlation radius. 

A comparison the of values of the exciton binding energy and the Coulomb interaction 
at distance d, e^/ nd, in Table IIIII leads to the conclusion that the correlations between 
excitons can be neglected only when a significant part of them is dissociated. However, 
in this case the concentration n that controls the blue shift in Eq. fll.ll) is not the exciton 
concentration but the sum of the exciton concentration and concentration of electrons or 
holes. 

The existence of region V is the result of the tail of the exciton - exciton repulsion 
potential. In case of a short range potential (e.g., hard circles) an overlap of the particle 
wave functions competes with the repulsion and the region of classical liquid behavior does 
not exist. This happens to excitons in one well where there is no dipole - dipole repulsion 
and an increase of the concentration leads to the Mott transition but not to a classical liquid. 
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FIG. 8: Different role of the exciton - exciton correlations in an exciton system. The shaded area 
shows the part of the region V where quantum correlations set in. 

TABLE III: Estimates for the parameters of Figl8]for GaAs/AlGaAs quantum well structures with 
different values of L. 



d (nm) 


(K) 


eV^d (K) 




(K) 


6V4(i4 (cm-2) 


l/d^ (cm-2) 


10 


54. 


139 


3.8 




2.2 X 10^° 


1.0 X 10^2 


12 


49. 


116 


1.8 




1.1 X 10^° 


0.69 X 10^2 


14 


44. 


99 


1.0 




0.58 X 10^° 


0.51 X 10^2 



Two of the lines in FiglHl between regions III and IV and between regions II and V, actu- 
ally comprise one line at which the gas parameter condition is violated. The gas parameter 
is the product of the concentration and the scattering crosssection squared. Between regions 
II and V this crosssection is classical, ~ tq, while between regions III and IV it is quantum, 
~ l/kx- The other two lines, between regions II and III and between regions IV and V, 
separate classical and quantum interactions between excitons. Also, it is necessary to note 
that at the bottom of region V some quantum coherent phenomena are possible. 

It is important to note that the lines separating different regions in FigJH] do not cor- 
respond to sharp transitions. Crossing of one of the lines by changing the temperature 
or the concentration leads to a gradual change of the correlation between excitons. FiglH] 
demonstrates only the role of correlation but not phases of the system. 
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VII. CONCLUSIONS 



We studied an exciton system in coupled quantum wells where electrons and holes are 
confined in different wells and the main interaction between excitons is a dipole - dipole 
repulsion. We found that in the most part of the temperature - concentration plane the 
system is characterized by a strong exciton ~ exciton correlation. At some parts of this plane 
the system behaves as a gas with a strong pair correlation. In other parts where the gas 
parameter condition is violated, i.e., where the probability to find more than two excitons 
close to each other becomes of the order of unity, the correlation is multi-excitonic and the 
system has to be considered as a liquid. In particular, at low concentrations degeneracy of the 
system is accompanied by a setting in of multi-exciton correlations. At high concentration 
the strong confinement of each exciton wave function due to repulsion between excitons 
suppresses quantum correlations. The blue shift of the exciton luminescence has a different 
value and a different dependence on the exciton temperature and concentration depending 
on how close excitons can come to each other. Therefore it is a sensitive tool for measuring 
of the exciton - exciton correlations. 
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APPENDIX A: EXCITON - EXCITON SCATTERING 

Here the problem of the exciton - exciton scattering is considered under the condition of 
small energy of relative motion, 

E = -TT- < —7 , Ala 
M Kd ^ ' 

and large wavelength, 

k < h/2(f . (Alb) 



For two excitons the center of mass momentum and the momentum of relative motion 
are defined as 

K = fei + fe2 , = . (A2) 
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The two-exciton wave function is factorized 

^{r^, r,) = 4= e^^(^^+^^)/2^(ri - r,) , (A3) 
V S 

and the Schrodinger equation for the wave function describing their relative motion is 

- — V^ij{r) + U{r)ij{r) = Ei){r) . (A4) 



Under the condition Eq. ( lAlal) the minimal distance between excitons is much larger than 



d (i.e., at r ~ c? the wave function is negligibly small) and the interaction potential can be 
approximated with Eq. (l3.7p . The condition of a very long exciton wavelength, Eq. (1Alb[) . 
makes it possible to simplify Eq. flA4p in two regions. In the region where the distance 
between excitons is much smaller than the wavelength of their relative motion, kr <^1, the 
coordinate dependence of the wave function comes only from the potential energy and the 
characteristic scale r ~ (P /h. At this scale the kinetic energy can be neglected and Eq. (]A4l) 
is reduced to 

VV(^) - i^{.r) = , A;r < 1 . (A5) 

When the distance between the excitons is much larger than tq = [e^d? / kEY^^ = 
{2d'^ /k'^by^^ the interaction energy is small compared to the kinetic energy and 

V^^/'(r) + k'^ip{r) = , r > tq . (A6) 

Due to krQ = (kd'^/bY^^ <^ 1 the two regions overlap at tq ^ r ^ l/k. 

Solutions to both Eqs. flASI) and flA6P are expressed in Bessel functions. When kvQ -C 1 
only S scattering is important and it is enough to find angular independent solution of 
Eq. (lA5p . The solution that goes to zero at r — > is 

ij{r) = A,Ko (^) = A,Ko (2kr, .{^\ , fcr « 1 . (AT) 



'br 

The solution to Eq. (lA6p describing scattering is 



^{r) = e^*'--^ + A2H'i\kr) . (A8) 



Making use of asymptotes [35| 



Kq{z) = -\xi^-C + 0{z^hiz) , 1^1 <1, (A9a) 
H^^\z) = l + -{\r,'± + c]+0{z^\nz) , k| < 1 , (A9b) 
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where C = 0.577 is the Euler constant it is easy to match the solutions in the intermediate 
region vq <^ r <^ 1 / k: 



1i / , kr ^ 
1 + - 1„- + C 



(AlO) 



This gives 



A, 



m/2 



(All) 



ln( A;d2 /h) + 3C - In 2 - Z7r/2 ' ln( A;d2 /h) + 3C - In 2 - i7r/2 

According to Eq.l lAlbl) the argument of the logarithm is small and with the logarithmic 

accuracy 



A, 



\Yl{k(P/h) ' 



A,, 



m/2 



lYiovd the asymptote [35| 
it follows 

7/;(r) = e**'^ + 




^i(z-K/i) 



in{kdyb) ■ 



\z\ > 1 



Akr 



fcr > 1 



ln(fcdV6) V 2fcr 

That is the scattering crosssection (in 2D case it has units of length) is 

^2 



TT 



cr 



k\n\k(P/h) 



(A12) 



(A13) 



(AM) 



(A15) 
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